THE ENERGY-SCALING APPROACH TO 
PHASE-ORDERING GROWTH LAWS 
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We present a simple, unified approach to determining the growth law for the characteristic length 
scale, L{t), in the phase ordering kinetics of a system quenched from a disordered phase to within an 
ordered phase. This approach, based on a scaling assumption for pair correlations, determines L{t) 
self-consistently for purely dissipative dynamics by computing the time-dependence of the energy 
in two ways. We derive growth laws for conserved and non-conserved 0{n) models, including two- 
dimensional XY models and systems with textures. We demonstrate that the growth laws for other 
systems, such as liquid-crystals and Potts models, are determined by the type of topological defect 
in the order parameter field that dominates the energy. We also obtain generalized Porod laws for 
systems with topological textures. 
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I. INTRODUCTION 

The quench of a system from a disordered phase to 
an ordered phase is a non-equilibrium process in which 
energy is dissipated and topological defects, if present, 
are eliminated. Typically, the system develops a scaling 
structure with a single length scale that evolves in time 
as the various broken-symmetry phases compete to select 
the ordered phase [Q-^. In the thermodynamic limit, this 
length scale will grow without bound, and the scaling 
structure will hold at late times. When scaling holds, 
any theoretical or experimental analysis of a system is 
simplified. It is then natural to explore the growth law 
of the single length scale, L{t). That is the aim of this 
paper, aspects of which have appeared elsewhere . 

Previous investigations of growth laws have been car- 
ried on a case-by-case basis, and many of the predictions 
have been controversial. The approach presented here 
is simple yet powerful, dealing with all phase ordering 
systems within a common framework. This approach is 
based on the role played by topological defects in fix- 
ing the large-momentum behavior of two-point correla- 
tion functions. The only restriction is that the dynam- 
ics be purely dissipative, the major assumption that the 
scaling hypothesis is valid. Even this assumption may 
be relaxed: in systems with more than one characteris- 
tic scale, our approach sets one relationship between the 
length scales. 

Systems with scalar order parameters, such as bi- 
nary alloys and Ising models, have been well studied 
p]-p|,p|-[20[ . In such systems domains of both phases 
grow, and the intervening domain walls (the characteris- 
tic topological defects of scalar systems) decrease in to- 
tal area and hence dissipate energy. Recently there has 
been a growing interest in systems with vector and more 
complex order parameters [ p|. A series of experiments 
[Eh-E3[ and simulations p4|-^f| have explored such sys- 



tems. A characteristic aspect of many of these systems is 
the existence of stable topological defects with singular 
cores. These defects include domain walls, vortex lines, 
and points [0. If a system has a scaling structure, then 
the single length-scale L{t) will also characterize the cur- 
vature of any domain walls and vortex lines, and the sep- 
aration of any point defects. One can use the nature of 
the defects to determine the short distance correlations in 
the system. This connection between the topological de- 
fects of the system and its correlations can be exploited. 
The scaling assumption and the knowledge of the defect 
structure are together sufficient to determine the growth 
law. 

We calculate the growth law of the characteristic 
length scale, if scaling exists, for quenched systems with 
either scalar or vector fields, with either short- or long- 
range attractive interactions, and with or without con- 
servation laws. We do this by considering the time- 
dependence of the energy density as the system relaxes 
towards a ground state. First we calculate the energy 
density of the system. Then we equate its time-derivative 
to the rate of energy-density dissipation independently 
calculated from the local evolution of the order parame- 
ter. From this we self-consistently determine the growth 
law, L{t). Our approach is independent of the details 
of the system and of the initial conditions, and hence 
reflects the observed universality of growth laws among 
physical systems and simulations. 

The symmetries of the system are reflected in the topo- 
logical defects seen in the system. The defect structure 
determines the asymptotic behavior of correlations, so 
that we do not need to use any approximation schemes in 
our approach. For short-range interactions we find three 
regimes (shown in figure |l]) : defect dominated scalar and 
XY [0(2)] systems both with and without relevant con- 
servation laws, and spin-wave dominated | ^2| systems in 
which conservation laws are always relevant. In the de- 
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feet dominated regime topological defects provide short- 
scale structure that determines the energy density and 
the growth laws. In the spin- wave dominated regime, 
topological defects, if present, do not dominate the en- 
ergetics or determine the growth law. For all regimes, 
we find power-law growth laws which are independent of 
the spatial dimension of the system, apart from special 
cases of one-dimensional (ID) scalar systems and two- 
dimensional (2D) XY models. For the marginal cases 
separating regimes we find novel logarithmic factors in 
the growth laws. 

It is informative to contrast our approach with 
renormalization-group (RG) work by Bray (see 
also ^,^). In common is the universality with respect 
to the details of the phase-ordering system, and the as- 
sumption of dynamical scaling. There are three signifi- 
cant differences. The first is that the RG approach only 
determines the dynamical exponent z, i.e. the power law 
1/z of the growth, but does not determine any logarith- 
mic factors. The second one is that the RG approach 
implicitly assumes the existence of a simple fixed point in 
which all correlation functions scale with a single length, 
while our approach only requires the scaling of the two- 
time two-point correlation function. The third difference 
is that the RG approach is limited to systems with a 
relevant conservation law, while our approach also ap- 
plies with non-conserved dynamics. In contrast to the 
RG methods, our "Energy-Scaling" approach places the 
growth laws of both conserved and non-conserved scaling 
systems on equivalent theoretical footing, and explicitly 
predicts any logarithmic factors. In principle, the results 
of this work could be reached by the RG approach us- 
ing an explicit RG structure for non-equilibrium phase- 
ordering problems, but, to our knowledge, that has not 
yet been developed. 

In section |l| we introduce the model. In section III we 
obtain growth laws for a variety of 0(n) systems. We 
introduce the scaling assumptions and then calculate the 
energy density and the energy-density dissipation. We 
also discuss systems with long-range attractive interac- 
tions. These results have mainly been presented before 
M, though in a brief form (see also pi). In section IV 



we discuss systems without topological defects, then con- 
sider the collapse of isolated defect structures such as 
spherical domains, 2D XY point defect pairs, and non- 
singular textures. We use these results to obtain growth 
laws for 2D XY models, and also asymptotic correlations 
and growth laws for systems with topological textures. 
These have mostly not been presented before and pro- 
vide a detailed perspective on scalar, XY, and textured 
phase-ordering systems. In section ^ we review relevant 
work on phase-ordering growth laws. We then discuss 
our results in section ^ and summarize them in section 
Vil- 



li. THE MODEL 

A generic energy functional for a phase-ordering sys- 
tem with an n-component order parameter, 0(x), and 
short-range interactions, is 



(1) 



where d is the spatial dimension and V{(j)) is an isotropic, 
"Mexican-hat" shaped potential such as 



(2) 



After a temperature quench into the ordered phase, the 
equation of motion for the ordering kinetics of the Fourier 
components 0k [|47| is 



9* 0k 



(3) 



where we only consider systems with purely dissipative 
dynamics. We work at temperature T = 0, with no ther- 
mal noise, since T is an "irrelevant variable" within the 
ordered phase |4^,^^. The conventional non-conserved 
model-A and conserved model-B dynamics are fi — 
and fi = 2, respectively. However, any /i > enforces a 
global conservation law for Q . 

Shortly after the quench, the magnitude of the order 
parameter nearly saturates and evolution takes place by 
the motion of topological defects and, for vector systems, 
by the relaxation of the director-field of the order pa- 
rameter. The role of the potential, V {(/)), in these later 
stages of evolution is different between scalar and vector 
systems, thou^ in neither case does it dominate the gra- 
dient term in (yj). Near defect cores, where the order pa- 
rameter vanishes in order to reduce the gradient energy, 
both the gradient and potential terms contribute to the 
energy-density in proportion to the defect-core density. 
Away from defect cores, the field is close to saturation 
and varies over scales of order L{t), the characteristic 
length. Consider a small region, away from defect cores, 
with a uniform gradient in the director field and a uni- 
form magnitude |0| — 1^6. The local energy density, 
H = (V0)^ -I- V{(j)), will be of the same form for any 
potential V{4i) with a quadratic minimum: 

H[0]~(1-<5)2(V0)2+4K,<52. (4) 

The local energy density is then minimized if 

S^{V^r/{Wo). (5) 

Since V0 ^ 1/L(t), the contribution to e from the po- 
tential term, of order 6'^ ~ is clearly subdominant 
for large L{t). Including non- uniform gradients in the di- 
rector field, and corresponding gradients of S, only con- 
tributes more subdominant terms. For vector systems, 
the energy due to the magnitude variations of the order 
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parameter balance between the potential and gradient 
terms, but both are dominated by the gradient term of 
the director field. For scalar systems, the director field 
is uniform away from defects and so the energy density 
is dominated by the domain wall energy, which is bal- 
anced between the potential and gradient terms. In both 
cases, the gradient term is proportional to the local en- 
ergy density. Thus the scaling behavior of the average 
energy density of the system, e = (7i), is captured by 



Lit) 



(10) 



(6) 



where is the momentum integral J d'^k/ {2ttY, and the 
angle-brackets represent an average over initial condi- 
tions, or, equivalently, a spatial average in the thermo- 
dynamic limit. We independently calculate the rate of 
energy density dissipation, e = dte, by integrating the 
contribution from each Fourier mode: 



(7) 



where we use the equation of motion (|^) to obtain the 
second line. We see the special role of the two-point 
equal-time correlation function. 



5(k, t) 



(8) 



which determines the energy density, and of the two-point 
time-derivative correlation function, 



(9) 



T(k,<) = (dtcj)k-M 



which determines the rate of energy-density dissipation. 
These correlation functions will be the basis of our dis- 
cussion of growth laws for systems that satisfy dynamical 
scaling in the following section and subsequently of iso- 
lated defects and growth laws in scalar, XY, and textured 
systems. We will equate the time-derivative of e from (||) 
to e from (^) , and self-consistently determine the growth- 
laws. 

Before proceeding, we digress to discuss what can be 
learned from dimensional analysis of the dynamics in Eq. 
(^. From the linear part, associated with the gradient 
term in H, of the equation of motion (|^) we have the di- 
mensional relation [L] — What other length 
scales are there in the problem? The nonlinear part 
of (H), associated with the potential term in H, gives 

[L] = [^] where ^, given by V^^^"^ for the potential (||), is 
the core size of any topological defects present, e.g. the 
width of a domain wall. Finally, there is the length scale 
^0 associated with any short-range correlations in the en- 
semble of initial conditions. Thus dimensional analysis 
gives [L] = [ti/(2+M)] = [^] = implying 



The growth law can only be changed from the naive 
(without ^ or ^o) dimensional result L ~ ^1/(2+/^) non- 
trivial behavior of F at late times, when its arguments 
approach zero. However, cannot enter into the growth 
law unless scaling is violated (see, e.g. | |25| , p8| , p9| ) . If 
did enter into the asymptotic growth law, then chang- 
ing the initial condition from the state at i = to one 
at ii > (with a larger ^o) will lead to a multiplicative 
change in the amplitude of the growth law, rather than 
merely a shift of origin of the time coordinate. This is un- 
physical, so any dependence implies a lack of scaling. 
Conversely, systems which break scaling need the initial 
conditions to set, through the initial correlation length, 
the relative amplitudes of the multiple growing length- 
scales. [The earlier argument against dependence does 
not apply to systems which break scaling, since changing 
^0 in such systems is no longer simply equivalent to an 
offset in the time-coordinate, but also changes the rela- 
tive amplitudes of the growing length scales.] It is pos- 
sible, but unlikely, that multiple length scales will have 
relative amplitudes set only through the core scale ^. So 
we expect scaling violations if and only if we observe 
dependence in the amplitude of the growth law. For scal- 
ing systems, therefore, we may set the second argument 
of the function F in ( pX| ) to zero. Then, in cases where 
no topological defects are present, or the defects do not 
dominate the dynamics, L{t) will not depend on ^ either, 
and ( p^ ) gives L{t) ^ ^1/(2+*') This is the generic result, 
covering the region n > 2 in figure ^. For n < 2, the 
defects do dominate the dynamics (in that they domi- 
nate the energy density, as we shall see in section III B| ) , 
and L{t) can depend on the core scale ^. The remainder 
of the paper is mainly devoted to understanding these 
cases. 



III. GROWTH LAWS 

A. Scaling assumption 

Many phase-ordering systems are empirically found to 
scale at late times from a quench, after initial transients 
have decayed. Accordingly, the correlation function of 

the order parameter, C(r, t) — ^0(x, t) ■ (j){:ic + r, t)^ , ex- 
hibits the scaling form C{r,t) = f{r/L{t)), with a single 
characteristic length scale, L{t), and a time-independent 
scaling function f{x). This is the dynamic scaling hy- 
pothesis | [Tl| . Fourier transforming the scaling form of 
the correlation function, we obtain the scaling form of 
the structure factor. 



S{k,t) = L{tfg{kL{t)), 



(11) 



where d is the spatial dimension and g(y) — 
J d'^xe'^^'y f{x). This scaling hypothesis can be gener- 
alized to two-time correlations by dimensional analysis. 
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= k-'^g{kL{t),kL{t'),t/t'), 



(12) 



so that g{y) = y '^g{y,y,l). From this we obtain the 
scahng form of T(k, t), 



dtdt' 



Mt) ■ 0-k(t') 

L^k^^^hikL), (13) 



where h{x) is a new scahng function | ^0| . 

Our Energy-Scahng approach is based on the scahng 
behavior in equations ( pi] ) and (|l^). Since the structure 
of the system determines the energy, we can obtain the 
energy as a function of L. The time dependence of the 
energy wiU come solely through £(t), since the structure 
is invariant up to rescaling. However, the rate of energy 
dissipation is directly determined by T(k, t) and can be 
calculated as a function of L and L, with no other time 
dependence. By equating the time derivative of e from (|^) 
with e from (|^), we self-consistently obtain i as a func- 
tion of L. From this we solve for the growth law of the 
length scale, L{t). We first consider systems where n < d, 
which have singular topological defects ^l[] . We will dis- 
cuss sy stems without topological defects (n > d + 1) in 
section EVA and systems with non- singular topological 
textures {n = d+1) in section IV C . 



B. Energy density 

To calculate the energy-density we use the scaling form 
of the structure function (pT|). If the thermodynamic 
limit of the phase-ordering system exists, then the infra- 
red (IR) limit of the integral (||) is well behaved. Hence 
either the momentum integral converges in the ultra- 
violet (UV) and e L^^ is extracted using the scaling 
form and a change of variables, or momenta on the order 
of the UV cutoff dominate the integral. 

When structure in the UV limit dominates the integral, 
we need to know the behavior of the correlation function 
in that limit, with kL 3> 1. This small-scale structure 
will only come from topological defects, since small-scale 
non-defect structures relax quickly via the dissipative dy- 
namics. As a result, the structure factor is proportional 
to the density of defect core, 5(k, i) oc pdof for kL ^ 1. 
Since in an n-component model in d-dimensions the de- 
fect core will have dimension d — n, it follows from scaling 
that the core density pdcf ^ L'^-'^/L'^ ~ L"" Hence 
5(k, t) ^ i^", and the scaling form ( pi] ) implies [ pT| - 



5(k, t) - L-" k 



-(d+n) 



kL > 1. 



(14) 



The form of the Porod tail is purely geometrical in ori- 

tand docs not depend on the overall defect structure 
, or on the details of the defect core (this holds even 



for defects with asymmetric time- independent core struc- 
ture, see 1^^). This generalized Porod's law is valid when 
n < d, so that singular topological defects exist. 

Using the asymptotic expression (^4|) as needed in the 
integral for the energy density (^), and imposing a UV 
cutoff at fc ^ 1/^, we obtain |5J] 



i-2 ln(L/b , 
L-' , 



n<2 , 
71 = 2 , 
71 > 2 . 



(15) 



The integral (^ is UV divergent for tt, < 2 and conver- 
gent for n > 2. For tt, = 2, when the integral is logarith- 
mically divergent, we impose a effective lower cutoff at 
k ^ 1/L, which is the length scale at which ( p^ breaks 
down and which reflects the relative lack of structure at 
scales longer than L{t). We see that the energy is dom- 
inated by the defect core density, pdof, for ?i < 2, by the 
defect field at all length scales for 77 = 2, and by varia- 
tions of the order parameter at scale L{t) for tt, > 2. 



C. Energy-density dissipation 

We calculate the scaling behavior of the rate of energy- 
density dissipation, e, as a function of L{t) in an analo- 
gous manner. If the energy-dissipation integral (|^) con- 
verges we use the scaling form ([l^ ) and change variables 
to obtain e ~ L'^~^L^ . For the cases when the momen- 
tum integral diverges in the UV, we must evaluate the 
time-derivative correlation function in the kL ^ 1 limit. 

In the UV limit, structure comes from defect cores. 
The field sufficiently close to a defect core comoves with 
the core. If the core has a local velocity v, then 



(16) 



close to the core (i.e. at distances small compared to L), 
where i?^ is a rotation matrix. In the momentum repre- 
sentation, the second term of ( p6| ) probes the structure of 
the defect and scales as kL, while the rotation term does 
not and only scales as L/L. Hence the rotation term is 
negligible when kL 3> 1. This gives dt4>\i kv4>k in that 
limit, and 



9*4 • 5t4k) - (v^), (4 • 4k) , 



2 — d— 71 j-^—r 



kL > 1, 
fcL>l, (17) 



where (t^'^)^ is the square velocity of defect core, aver- 
aged over the core elements that contribute to structure 
at momentum k. To obtain the second line we have used 
the generalized Porod's law ( p^ and have implicitly as- 
sumed that all lengths are large with respect the core 
size: k~^,L ^ ^. We emphasize that ( p^ is valid suffi- 
ciently close to a slowly moving defect core. Since dt4> is 
purely dissipative from (^, the defect core in (|l^) must 
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be moving "downhill" in the energy landscape of the lo- 
cal order parameter field and, for extended defects, under 
the influence of the surface or line-tensions of the defect 
core itself. 

The characteristic value of the defect velocity v is L. 
If the energy-dissipation is dominated by the evolution 
of large defects at the characteristic size L, then using 
Eq. (|l7| ) and replacing Iv'^) f, with within the energy- 
dissipation integral (0) is appropriate. We show below 
that this is the case for n < d or for n > 2. If the 
integral converges in the UV this is equivalent to using 
the scaling form ( p^ . If the integral diverges in the UV, 
this implies that the energy dissipation is due to the slow 
evolution at characteristic scales (such as shrinking do- 
mains or vortex loops), leading to a reduction of the core 
volume. The actual asymptotics of ([l^) is in general not 
give n by u sing (i'^)^, i^, as we will see explicitly in sec- 



tion IV B, so we must directly check the appropriateness 
of the substitution in calculating e. 

For n < 2, when the energetics is dominated by the 
defects (following Eq. (^5|)), the rate of energy dissipa- 
tion can be written as the rate of change of the energy of 
all the defect features Q : 

d 

e^-y dln{l,t)E{l), 



jiom) 



dij{i,t) 



dE{l) 
dl 



(18) 



where n(L t) is the number density of defect features of 
scale / [pM, E{1) is the average energy of a defect feature 
of scale Z, and is the number flux of defect features. 
We use the continuity equation, dn/dt + dj /dl — 0, to 
obtain the second line in (|lj) . The total number density 
of defect features, N , scales as the inverse scale volume, 
N ^ ^jE'^, and hence iV is slowly varying for times of 
order Ljh. Since defects only vanish at the core scale ^, 
we have N ~ i(C)- This implies that j(/) is independent 
of I for Z ^ L, in order to provide a steady rate of defect 
extinction. We would like to know whether the integral 
in Eq. ( p^ ) is dominated by small scales, with I L, 
and so we need to know E{1) in that limit. We assume 
here that small defect features have an average energy 
related to the scahng form (na) for the energy-density: 



E{1) 



n<2, 

/''-^ n > 2. 



(19) 



This is ex plicitl y confirmed for isolated symmetric defects 
in section IV B , Using (|l^) , we see that for d > n or n > 2 
the integral in (18) is well behaved at I L and the 
integral dominates the JiO^iO term. For these cases, 
structures with scales and separations I ^ L{t) always 
dominate the energy dissipation, so that {v'^) j. ^ can 



be used in evaluating e with equation (|^. For d = n < 2 
the integral in ( |l8| ) diverges for I <^ L, and dissipation 
is dominated by defect pairs annihilating. We treat the 
case d — n — 1 below. For d = n = 2 the integral is 
logarithmically divergent, and e has contributions from 
defect pairs at all scales. We treat this case in section 



IV B by determining {v )k directly. 



For systems with 7i < d or n > 2 we use {v'^)^ 
in (|l7|), within the energy-dissipation integral (0). This 
gives 



L-"ln(L/Oi^, n- 



Ai < 2 
/i = 2 
/i > 2 



(20) 



apart from n ~ d < 2. For n + fi < 2 the integral is 
UV divergent and we have used the fcL 3> 1 form ( [TtI ) 
between its lower limit of applicability, k ^ L~^, and 
the UV cutoff at k ~ 1/^. For n + fj, = 2 the integral 
is logarithmically divergent. Again we use the kL 3> 1 
asymptotics in its region of applicability. [ The scaling 
contribution from kL ^ 1 merely changes the effective ^ 
in the logarithm of (20).] For n -f /i > 2 the integral is 
convergent and we use the scaling form (|l^). Equating e 
from ( po| ) to the time derivative of e from (p_5|) , we obtain 
the growth laws shown in figure |^. 

As an example, consider nonconserved scalar fields. 
For this case, with /i = and n = 1, Eq. dl3) gives 
e ~ implying e ^ — L/L^^, while ( |20| ) gives 

e — i^/L^. Equating these two results for e yields 
L ^ 1/L, implying L ~ t^^^. The other results displayed 
in figure |^ were obtained in the same way. 

For ID scalar systems {n = d = 1), equal-time correla- 
tions scale but r(k, t) breaks scaling due to rapidly 
annihilating domain wall pairs [^. We can still apply 
the Energy-Scaling argument to the exponentially sup- 
pressed interaction energy between domain wall pairs, ex- 
pressed as an energy density, eint ~ e~^/^/L. We equate 
the time-derivative of this to the energy-dissipation from 
equation (|2^) (since the interaction energy changes from 
the slow (w^) ~ P evolution of the defects). This leads 
to logarithmic growth L(t) Int for both non-conserved 
and conserved scalar models in one-dimension, as shown 



in table [V, 

The growth laws shown in figure and table |l^ are 
only the leading time-dependence of the properly inte- 
grated equations of motion for the length scale. Further 
corrections come from subdominant parts of the energy 
integrals, and from subdominant corrections to the scal- 
ing forms of the correlation functions. Corrections to in- 
tegrals dominated by the UV limit are given by contribu- 
tions from kL 1, and by corrections to the generalized 
Porod's law (p^. The latter are not generally known, ex- 
cept for three-dimensional scalar systems where the lead- 
ing correction to Porod's law is |57|. Corrections 
to convergent or logarithmically divergent integrals can 
come from the UV cutoff. For instance, for non-conserved 
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XY systems the In L factors in (y_5|) and (^0|) will in gen- 
eral have different effective cutoffs, of order the core size. 
This leads to an additive correction to the growth law, 
0{L{t) / Xwt). It is important to be aware of such correc- 
tions for each particular system. 



D. Long-range attractive interactions 

For systems with long-range attractive interactions 
|p|, P^|j2^]4^ the rate of energy-density dissipation will still 
be given by (0), but the energy density has a new contri- 
bution 



(21) 



with < fj < 2. This reduces to the short-range case for 
(7 = 2, and is subdominant to the short-range interactions 
for cr > 2. This interaction can be motivated by a term 
in the energy-functional: 



LR 



(x -I- r) - (j){yi) 



(22) 



where we take tr > for a well defined thermodynamic 
limit. As in the short-range case, the local potential V{(j)) 
does not dominate the gradient, and can be neglected. 
We see this for < cr < 2 by generalizing equation (|^) to 
balance a long-range energy (1 — 5)^L~'^ with the poten- 
tial term AA5^ . The local energy density is minimized if 
i5 ^ l/L'^, where |^| = 1 — 6. This corresponds to a 
tail in the profile of an isolated defect, a generalization 
of the l/r^ tail for vector systems with only short-range 
interactions. This tail does not change the generalized 
Porod's law (|lj) for either short- or long-range interac- 
tions, because the local structure of the defect does not 
change as the length scale L{t) grows. As before, the 
structure factor for kL 3> 1 is simply proportional to the 
density of defect core. Using Porod's law in we find 




71 < cr, 
n = a, 
n > a. 



(23) 



Comparing with the energy-density from the short-range 
interactions (p^, we see that elr always describes the 
scaling of the energy density at late times for < a < 2. 

The asymptotics of T(k, t) are still given by (p7|). Ap- 
plying Eq. (18) to the long-range case shows that we can 
use {v'^}f, ~ within the energy-dissipation integral (|^) 
for all cases except rt = d < cr [Q. Comparing the rate 
of energy dissipation, still given by (pO|), but excluding 
n = d < a, and the time derivative of (p3|), we find L 
and hence L{t). The results are summarized in figure]^. 
We treat the cases n — d=l<a<2, where (w^)^, is 
dominated by small defect features, elsewhere and find 
L(i) ^ g. 



IV. SPECIAL CASES 

A. Systems without Topological Defects 

Systems with no topological defects [such as 0{n) mod- 
els with n > d -|- 1] will have no power-law UV structure, 
and hence will have convergent energy integrals, giving 
L{t) ^ t^/i'^+t^) ^ where a = 2 for short-range interactions. 
This result for L{t) also follows directly from the dimen- 
sional analysis discussed at the end of section ||, if the 
system scales. Any power-law UV structure would indi- 
cate a scaling density of singular structure — which we 
do not expect if locally stable defects are absent. For in- 
stance, Rao and Chakrabarti find stretched exponential 
tails to the asymptotic structure factor for 2D 7i = 4 
systems. Of course, this result applies to any system 
with initial conditions that evolve into a scaling structure 
without topological defects, and not just to ri > d + 1. 
For example, if the 2D XY model is quenched from be- 
low Tkt then the bound vortex pairs quickly annihilate 
and the scaling state has no topological defects (for a 
quench to T = 0) [pSj. In the non-conserved case scal- 



,1/2 



as 



ing is obeyed and the growth law is L 
expected. In the conserved 2D XY model quenched from 
below Tkt, t^^^ growth (or, more generally, will 
obtain, provided scaling holds. Of course, scaling can be 
violated even without topological defects — as seen in 
the conserved spherical model [p9|. 



B. Isolated Defects 

In this section we obtain the size l{t), as a function of 
time to collapse t, of an isolated symmetric defect. We 
will then use this, for scalar and 2D XY systems, to ob- 
tain the proper asymptotics of T(fc,t), i.e. an expression 
for (f^)fc in (|l7|). For scalar systems this will confirm the 
growth laws we have obtained. For 2D XY systems, we 
will be able to obtain the scaling growth laws (see ta- 
ble [V ) where previously we were prevented by our crude 
knowledge of (u^ ) ^ . [In the following section we will apply 
a similar approach to systems with non-singular topolog- 
ical textures.] 

Applying equations (||) and (^ to isolated defects, in 
a "microscopic" approach, is appealing but not trivial. 
For scaling systems, with many defects, the disordered 
initial conditions will cut off any momentum integral in 
the infra-red (IR). Equivalently, there is no structure at 
scales much greater than the inverse of the characteris- 
tic length, , which thus provides a natural IR cutoff. 
On the other hand, for solitary defects (such as a sin- 
gle spherical domain), an IR cutoff is provided by the 
structure and evolution of the single defect itself. We 
will apply the quasi-adiabatic assumption to the defect 
structure: that the order-parameter field is locally equi- 
librated to the defect core configuration. The breakdown 
of this assumption away from the core, due to a finite core 
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velocity, will provide a natural IR cutoff — cutting off the 
structure at small fc, corresponding to large spatial scales. 
We must apply the largest appropriate IR cutoff, and we 
will focus on small defects within larger phase-ordering 
systems, where both cutoffs to the IR behavior are appro- 
priate. Whenever possible, we check our approach with 
established results. We restrict our attention to the col- 
lapse of a spherical scalar domain of radius I with d > I, 
and the annihilation of a point defect/anti-defect pair of 
separation I in a 2D XY model {n = d = 2). We calculate 
both the structure factor S'(k, I) and the time-derivative 
correlations r(k, I), use these in the independent energy 
equations (H) and (j^), and self-consistently solve for 

We first consider generic systems with an isolated, 
symmetric, singular topological defect structure {n < d) 
Plf characterized by a single length-scale much larger 
than the core size, I ^ ^. For example, spherical do- 
main walls, or a point defect/anti-defect pair. The de- 
fect size, I, corresponds to the local radius of curvature 
for n < cZ or to the separation of point-defects for n = d. 
The structure-factor at small scales, kl ^ 1, is propor- 
tional to the core volume of the defect, l"^"", and has 
the asymptotic momentum dependence fc- (''+") of a 
"flat" or stable defect structure. [This asymptotic mo- 
mentum dependence leads to a generalized Porod's law 
for scaling systems [^-Q.] For kl ^ 1 the field around 
an element of defect core will comove with it [without 
a rotation term because of the symmetry of the defect] : 

dt4> = {dl/dt){V if)) radied, whcrc we have taken the radial 
component of the gradient. We then have, for the iso- 
lated defect, 

r(k, /) = {dl/dt)^ f-"" fc-(''+«-2) ^ (24) 

where the scaling functions (x) and fri^) tend to con- 
stants for large x. From (g), the energy of the isolated 
defect of scale I :s> ^ is given by 



1 

d''kl''-"k-^'^+"-^^ fs{kl) , 



(25) 



where we cut off the integral above at the inverse core 
size. The energy integral ( |25[ ) always converges at small 
k, as will be checked on a case by case basis. Evaluating 
the integral gives 



n = 1, 
n = 2, 
n>2. 



(26) 



[We include n > 2 for completeness. Isolated defects 
with n > 2 will generically be asymmetric [55 6^, and 
not described by a single scale I] The results for n = 1, 2 
are, of course, familiar: for n = 1 the energy is just 



proportional to the surface area of the domain, while 
the factor is an estimate of the surface tension; for 
n = 2, the energy is given by the defect core volume 
multiplied by a logarithmic factor from the 'far field', 
with a short-distance cut-off at the core scale ^ and a 
large-distance cut-off at the size I of the defect, reflecting 
the self-screening of the defect (a vortex loop for d — 
or vortex-antivortex pair for d = 2) at large distances. 

Next we evaluate the analogue of equation (|^), the 
energy dissipation rate: 

E{1) = - / d'^kk-^TikJ) 
Jo 

d'*fcfc-T(k,0- / rf'^fc A;-^ r(k, 

(27) 



= EiRil)+Euv{l) , 



where Ejji and Ejjv correspond to the momentum inte- 
grals over k < and k > respectively. The kl ^ I 
behavior, given by equation ( p^ ) with /^(fcZ) = const., 
can be used to calculate Ejjv- However, we will see 
that, unlike in the calculation of E{1), the IR behav- 
ior (through Em) sometimes dominates E{1), and so we 
need to determine both the IR structure of T(k, /) and 
the effective IR cutoff, kcut- For scalar systems, Ejr is 
only signiflcant for conservation laws with ^ > d, and 
fccut depends on whether or not the isolated domain is in 
a larger phase-ordering system or not. For 2D XY sys- 
tems, is always significant, and fccut is provided by 
the breakdown of the quasi-adiabatic approximation at 
large distances due to the non-zero evolution velocity of 
the defect core. 



1. Scalar Systems 

In scalar systems with d > 1, we consider a spheri- 
cal domain with radius I. Under the quasi-adiabatic ap- 
proximation, the field inside and outside the domain are 
taken to minimize the energy subject to the configura- 
tion of the domain wall. In the IR limit, with kl ^ \, 
the field is dominated by the saturated value inside and 
outside the domain walls and so ~ I'''. For conserved 
fields we know that 4'k=o does not evolve with / and so 
the quasi-adiabatic approximation must break down as 
fc — > 0. The finite evolution rate of the domain will de- 
termine the effective kcut- Using 0^ ~ Z'', we easily obtain 
the IR asymptotics: 



5(k,/) 

r(k,o ~ {di/dt)'^P'^-'^, 



kcut < fc < / \ 
kcut < fc < r\ (28) 



where we have not yet determined kcut- This IR limit, 
combined with the UV behavior from ( p4[ ) and the appro- 
priate kcut, will be sufficient to determine the evolution 
of an isolated domain. Note that the small kl forms (pa) 
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are consistent with the general scahng forms (p4[). In 
principle the complete scaling functions, fs and fx, can 
be calculated. In practice, however, we do not need the 
precise structure for /cZ 1, as its scaling properties with 
I are embodied in the IR and UV limits. 

For non-conserved systems, the IR asymptotics imply 
that the energy-dissipation integral ( p7| ) converges in the 
IR and is dominated by Ejjv , leading to 



E ■ 



'{dl/dtfl'^'^C 



(29) 



Comparing this with the time derivative of ( p6|) , recovers 
the standard result j|] dl/dt ^ —l^^ for d > 1, which 
gives 



1/2 



(30) 



where t = ^max — t is the time left before the defect an- 
nihilates. 

We can also extract a scale-dependent effective mobil- 
ity for non-conserved systems p2[, r]{l), by taking 



dl_ 



dE/dl 



(31) 



so that ri{l) is the proportionality factor between the driv- 
ing force per core volume and the evolution speed. This 
gives rj £^ for n = 1 < d. While the mobility of a scalar 
system depends on the details of the potential through 
^, i.e. on the surface tension, the evolution velocity dl/dt 
is independent of those details (for free energy function- 
al of the form (|]) ) ; this was first observed by Allen and 
Cahn II . 

For conserved systems with fi < d, the IR limit still 
does not dominate E and the dependence on I is still 
captured by Euv ■ The results are summarized in table |: 
for /i < 1 the UV limit dominates the energy-dissipation 
integral; for /i = 1 the integral is logarithmically diver- 
gent and is cut off by the detailed structure at fcZ ~ 1; for 
fi > 1 the integral of E converges in the UV, and struc- 
ture at kl ^ 1 dominates. The scaling behavior of this 
contribution is given by power counting on the second 
Une of (H). 

For strong-enough conservation laws, fJ. > d, the 
energy-dissipation integral no longer converges at small 
k. It is dominated by the IR behavior of Ejn and we need 
to know the appropriate cutoff /ccut • First we consider the 
breakdown of the quasi-adiabatic approximation far from 
the defect. This provides the only cutoff for a solitary do- 
main in an infinite system. The cutoff is provided by the 
behavior of the field (f> far from the defect, where we lin- 
earize the evolution equation around the ordered state 
= 1 and can neglect all but the lowest derivatives. 
From equation (||) we have 



kl < 1. 



(32) 



Any breakdown of quasi-adiabaticity (i.e. the field not 
following the static configuration given by the configu- 
ration of the domain wall) is due to the finite evolu- 
tion speed dl/dt. Dimensionally from Eq. (^, dl/dt ^ 



Ak^^y ^ , which determines a velocity-dependent IR cutoff 
/cvoi ^ {dl/dtY/'^'^~^\ This cutoff can be used for kcut in 
Eq. (g^ to determine the evolution of an isolated domain 
(with I ^ £/) in an infinite system. However, within the 
context of a larger phase-ordering system, other evolving 
domains provide a natural cutoff /ccut ^ L^^ , where L is 
the characteristic length-scale of the system. This cut- 
off can be used to determine the ^-dependence of dl/dt 
(i.e. the "scaling" behavior of the domain), again with 
Eq. (p7|). The two results can be compared and we 
find that the velocity cutoff is only larger than for 
domains with I < L^"^/"*. This is a vanishing propor- 
tion of domains in the scaling limit, so we have used 
^cut ^ in table ||. [The difference between the cut- 
offs does imply that for scalar systems with > c? a 
solitary domain will evolve qualitatively differently than 
a small domain in a phase-ordering system.] Comparing 
the dominant energy-dissipation {Em in this case) with 
the time-derivative of (pq), we find table Ifor d > 1. We 
obtain the standard I ^P^/^ for /i = 2 < c? |^ , and recover 
the ^ — d = 2 result of Rogers and Desai for small 
domains: dl/dt ~ —l/[P\n{L/l)], where, unlike ref. 
we have retained the defect scale in the logarithm. 

We can estimate (w^)j, for scalar systems with our re- 
sults in Table || This will enable us to determine the UV 
asymptotics of T{k,t) directly through ([l7|), and with 
that to check our previous results and the consistency 
of our approach. We assume that small defect features, 
with I <^ L, evolve as described in Table | with the lo- 
cal rms velocity of the defect core given by v{l) ~ dl/dt. 
We further assume that small defect features will evolve, 
at least in their I dependence, like segments of spherical 
domains, with a local radius of curvature decreasing as 
dl/dt, so that the number flux is given by j(l) = n(l)v(l). 
Then 



d—n 



^l^^^Jlv^{l)n{l)l 
Jdln{l)ld-^ ' 



kL > 1. (33) 



If small features (/ ^ L) dominate we use j{l) ^ ~ 
—LL"'^'^'^^^ in (p3|), while if small features do not domi- 
nate then the same approximation will give the expected 
(w^)fc ~ i^. To evaluate the integral, with ri = 1 for 
scalar systems, we take v{l) ^ dl/dt for / ^ L and use 
table H Only features large enough to contribute to the 
asymptotic structure of T(fc, t) are counted in {v'^)k, and 
this sets Irnin in (IH). Comparing Eq. (|l7| ) and j^) we 
see that /min = k~^ . We use this in equation (|33|) and 
summarize our results in table III. These results, used 
in equation ([l7|), give the true kL ^ I hmit of T(k, i). 
These asymptotics are not given by simply approximat- 
ing {v'^)k ~ L'^, at least for fi > d. However, the true 
asymptotics, used in (^, give the same growth laws as 
obtained earlier with our less detailed approach in sec- 
tion III (which has the advantage of requiring fewer as- 
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sumptions about the dynamics of small defect features 
within the phase-ordering system). For fj. < 1, where e 
is dominated by the UV asymptotics of T{k, t), then the 
proper asymptotics show that (w^)fe L^. Conversely, 
for ^ > d, where {v'^)k is not given by L^, we confirm 
that e is not controlled by the UV asymptotics of T(fc, t). 
Hence, the details of our treatment in section II], as well 



as the growth laws, are consistent with this more micro- 
scopic treatment. 



2. XY Systems 

We can apply the same quasi-adiabatic assumption to 
the 2D XY model: that the order-parameter field is given 
by the static configuration with boundary conditions im- 
posed by the defect configuration. For the 2D XY model, 
the hard-spin (Vb — > oo) energy is optimized by XY 
phases which satisfy V^O = away from the vortices. 
For an isolated defect/anti-defect pair separated by I, 
this leads to 



e{v) = tan" 



y 



x-l/2 



tan 



y 



1/2 



(34) 



where r = (x, y) and the order-parameter 0(r) — 
(cos 0, sin 0) p3. For quasi-adiabatic evolution, the time 
dependence is only through and the IR behavior of 
the spin-spin correlations is straight-forward to deter- 
mine from (p4): 



s'(k,^)--/Vfc^ fc;<ci, 

T(kJ) {dl/dt)^/k^, 



(35) 



Note once more that these results are consistent with 
the general scaling forms (^. Using (p5|), the energy 
integral converges in the IR and is given by (2g). For 
non-conserved systems, both Ejj^ and Ejjv (|27) are log- 
arithmically divergent. For conserved systems Em dom- 
inates and is infra-red divergent. For both cases we must 
determine the appropriate IR cutoff k^ut- 

The field far from the vortices can only can only follow 
their motion for infinitesimal dl/dt. For non-zero evo- 
lution speeds, the failure of the quasi-adiabatic approxi- 
mation far from the vortices provides an infra-red cutoff. 
Far from the vortices, the component of dtcf) parallel to 
the field is still dominated by Eq. (|2|), but the trans- 
verse component is given by the gradient term: {dt(j))j_ ~ 
— (— V^)''/^V^(^. Dimensionally from this equation, the 
transverse component provides a larger IR cutoff than 
Eq. (H), so that dl/dt ~ This breakdown of quasi- 

adiabaticity implies an IR cutoff fccut ^ {dl/dty^^^^^\ 
This cutoff can be used with the IR asymptotics ( |35| ) 
in the energy-dissipation integrals of Eq. (^) to deter- 
mine the evolution of a solitary vortex pair. We sum- 
marize the results in table ||. For closely separated vor- 
tex pairs within a phase-ordering system, other defect 



pairs will provide an effective cutoff at L~^. The largest 
cutoff will apply in Eq. (^^. We find that for non- 
conserved systems, the breakdown of quasi-adiabaticity 
provides the dominant cutoff for I ^ L/ln{L/^), and for 
conserved systems for I ^ L. Since in the non-conserved 
case fccut only enters logarithmically, for conserved and 
non-conserved systems any defect with I <^ L will have 
an effective IR cutoff provided by the breakdown of quasi- 
adiabaticity. For non-conserved systems, since the cutoff 
only enters logarithmically, the velocity screening is ef- 
fectively the same as the length-scale screening used by 
Yurke et al H], though from (H) we see that r(k, I) for 
the defect-antidefect is not screened in the IR through the 
static configuration. The effective mobility from equation 
^) is 7] - [In \£,dl/dt\]~^ for a small defect pair, which 
agrees with Ryskin and Kremenetsky ||6^. We see that 
for non-zero velocities the mobility does not depend on 
the system size, if the system is large enough. 

We can estimate (w^)^, for 2D XY systems with our 
results in Table ^ This will enable us to determine the 
UV asymptotics oiT{k, t) directly through ([T7|), and with 
that to determine the growth laws of 2D XY systems 
which scale — which we were unable to do in section III . 



We assume that vortex/anti- vortex pairs with separation 
I <C L evolve as described in Table |l|, with an rms speed 
of annihilation given by w dl/dt. Then the number flux 
of defect pairs is j{l) — n{l)v{l). If small features domi- 
nate then for I ^ L we use j{l) ^ j(^) ^ — LL^(''+^' in 
(^), using rt = 2, while if small features do not domi- 
nate then the same approximation will give the appro- 
priate {v'^)k ^ L^- To evaluate the integral we take 
v{l) ^ dl/dt for ; < L and use table 0. For 2D XY 
systems, all vortex separations contribute to the energy 
dissipation since T{\s.,l) has the same form for large and 
small kl. However we require the quasi-adiabatic ap- 
proximation to apply for the smallest separation, so we 
need k > kcut{lmin)- This leads to ^min ~ . We 
use Zniin in equation (^3|) and summarize our results in 
table III. These results can be used in equation ( p^ ) 
to get the true kL ^ 1 limit of T{\s.,t). We use these 
true asymptotics of T(k, t) in the energy-dissipation in- 
tegral (0). The energy-dissipation integral is logarith- 
mically divergent with e — LL~^ In (L/^). Equating 
this to the time derivative of e from ( [l5| ) we find that, 
although scaling is consistent, the growth law is not de- 
termined. However, imposing the scaling form (|l^) on 
the asymptotics ( |l7|) of r(k,i), using table IH, deter- 
mines the growth law, L ~ (i/lnt)^/^, for non-conserved 
systems and L ~ for conserved systems. These 

growth laws differ by logarithmic factors in comparison 
to XY models in d > 2 (see Fig. |l|), so that the 2D XY 
model is a special case of our both the treatment and the 
growth laws. 
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C. Textures 

Systems with n > d cannot have stable topological de- 
fects with singular cores, but systems with n = d+l can 
support a non-singular topological texture . We con- 
struct an isolated texture of scale X by stereographically 
projecting the field configuration from a c?-sphere of ra- 
dius R = X/2 surrounding an n-dimensional point defect 
in n-dimensional space onto the d-dimensional space of 
the physical system. For d = 1 we rest a circle (1-sphere) 
on the ID system. For d = 2 we rest a sphere (2-sphere) 
on the planar system. We set ^(x) parallel to the radius 
vector of the sphere at the point of intersection with the 
line joining x to the top of the sphere as shown 

in Fig. ^. This constructs an ri-component texture that 
winds once over the d-dimensional system. The texture 
is topologically stable and can only vanish when X ^ 0. 
We obtain 



(x) = (pxi,nx), 



(36) 



where the d-sphere touches the system at a; = 0. From 
Fig. m we have 



and 



Px — sin 2a, 

= 2xX/{X^ + x^), 



TT-x = — cos 2q!, 



(37) 



(38) 



The scale X of the non-singular texture is determined by 
nx = 0. We find 



(x) • ^(x + r) = 1 - 



2r'X^ 



(X2 + .t2) [X^ + (x + r)2 



■, (39) 



and 



S(k, X) = pdcf / d'^x / d''re*'^ "-(/.(x) • (/.(x + r), 



(2^)'^^(k) + 2pdcf^'v2 



2 ^ TT ^ Pdcf 



-2kX 



/cX>l, (40) 



where pdcf is the number density of defects (in the case 
of a single texture it is the inverse volume of the system) 
and we use the asymptotics of 



^{2^f/^X^-\kXf-'"^Ka,2-i{kX), 

(41) 



where K^(x) is a Bessel function of imaginary argument 
p5[ . For d = 2 this construction generates the mini- 
mal energy texture pa and has the same structure as 



unpublished results of Bray and Puri |Q . In general di- 
mension, this stereoscopic projection does not give the 
minimal energy, and we expect that a different convex 
surface (other than the d-sphere) is needed. However, we 
expect that the form of the structure factor for kX ^ 1 
will be unchanged and that the scaling form for a single 
texture will be 



5(k,x) = x2V(fc^) 



(42) 



which follows directly from the second line of equation 
(pO|). The time-derivative structure for kX ^ 1 is, using 
Eqs. (H) to dss"' 



T{k,X)=pdelX 



dXidX2 
d'^x f d'^r 



(x,Xi) 



•,d+2d+ 



-2kX 



kX > 1. 



This also leads to a scaling form for the time-derivative 
structure of a single texture: 



T(k, X) = X^X^'^-^hrikX) 



(44) 



We use these scaling forms to calculate the energy of 
an isolated texture of scale X. The energy will be, using 
equations (||) and (^^, 



E{X) = V 



d'^k 



/c^S'(k,X), 



[2ttY 

/ d'^kk^X^'^grikX), 
Jq 

^d^2 



(45) 



where V is the system volume, so Vpdd ^ 1 for an iso- 
lated texture. The rate of energy-dissipation is, using 
equations (||) and (|^, 



EiX) 



-V 



d'^k 



i2n) 



lT{k,X), 



d'^kX^X^'^-^hi 



(kX), 



_j.2^d-2_ 



Comparing these results, we see that 
X {d-2)X-^. 



(46) 



(47) 



For d > 2 an isolated texture will shrink (in agreement 
with Derrick's theorem |^8|) and vanish with X(t) ~ t^/^ 
as a function of time to annihilation, in d = 1 isolated tex- 
tures expand, while 2D textures are stable in this treat- 
ment. Indeed, in 2D Belavin and Polyakov present 
an exact solution for a static system of textures, with 
— 1 everywhere and an energy independent of the 
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scale and position of the individual textures. However, 
systems with mixtures of textures and anti-textures will 
be unstable |69| ]. 

We can use these results to discuss the phase-ordering 
of systems with many topological textures. However we 
can go far with the assumption that systems with tex- 
tures do not have more singular correlations than systems 
with point defects (n = d). With this, we obtain from 
(|l|) and (0) that S{k, t) < k"^'^, and r(k, t)< for 
kL ^ 1. We find that the integrals and ^ converge 
for systems with textures if d > 2. Hence for these tex- 
tured systems L ^ t^/C^+t^) ^ with or without conservation 
laws. This is in accord with the dimensional argument of 
section |l[ 

For non-conserved systems, we can use the results for 
isolated textures to calculate the asymptotic kL 3> 1 
structure factor due to a distribution of textures at small 
scales. Very small textures, with X <^ L, will collapse 
independently for d > 2. The flux of annihilating tex- 
tures is j(^) = ~N, and the number density of textures 
scales as an inverse volume N ~ L~'^. The flux of an- 
nihilating textures at small scales, is j{X) — nt{X)X, 
where nt{X) is the number density of textures at scale 
X. Using X ~ from @ and utiX) - N/X for 

X <^ L, we obtain 

S{k,t) ~ / dXnt{X)S{k,X), 
Jo 

^ j^j^-i2d+2)^ /tL>l, d>2, (48) 

and 

POO 

T{k,t)o^ / dXnt{X)T{k,X), 
Jo 

^^p-2d^ A:L>1, d>2, (49) 

where we have used the scaling structure of an isolated 
texture from equations ( ^2[ ) and (^). These results sat- 
isfy the bounds we obtained from the assumption that 
the structure will not be more singular than for systems 
with point defects — though the time-derivative corre- 
lations saturate the bounds. The k~^ tail predicted for 
the structure factor of three-dimensional n = 4 systems 
is in rough agreement with Toyoki's simulation results: 
S{k) ~ fc~7.5±o.2 ^j^g same system 

For d = 1 isolated textures will expand, because the 
energy of an isolated texture of scale X decreases with 
X ([47|). Hence we do not expect a singular annihila- 
tion process. This will also be true in d = 1 systems 
with conservation laws. Rather, we expect a non-singular 
combination of winding and anti-winding textures. The 
minimal texture scale will increase with time, and the 
exponential factors in ( pO[ ) and (^3|) will cause the mo- 
mentum integrals to converge in the UV. This results in 
a scaling prediction oi L ^ ^1/(2+^*)^ In fact, it can be 
shown that both non-conserved and conserved ID XY 
systems do not scale The scaling violations are as- 
sociated with the fact that the correlation length for 



the initial conditions enters in a nontrivial way. Extract- 
ing a length-scale from the energy-density, e ^ L^^, the 
non-conserved model has L ~ t^/"* while the con- 
served model has L ^ t^/^ [^-^. These growth laws 
are different than the scaling results, which is consistent 
with the lack of scaling. 

For d = 2, textures will either slowly shrink due to 
an instability caused by a soft-potential ||6^ , [7l|| , or will 
unwind with anti-textures in a non-singular annihilation 
process. In the former case, which does not seem to be 
the dominant process [ |69| , a similar approach to Eqs. 
( ^ ) and (|4^) finds bounds which lead to convergent en- 
ergy and energy-dissipation integrals for any annihilation 
rate. In the latter case, with non-singular processes, we 
also expect convergent integrals. In either case we expect 
L when scaling is obeyed — however, different 

growth laws and scaling violations are observed in sim- 
ulations with non-conserved dynamics [ ]69[ . Again, the 
scaling violations seem to be associated with a nontrivial 
role of the length scale characterizing the initial con- 
ditions 1^^. Conserved 2D n = 3 systems have not yet 
been investigated. 

The scaling assumption does not hold in the ID and 
2D systems with textures that have been investigated. 
Does it hold for d > 21 The scaling violations observed 
seem to be related to the weak interaction of the textures, 
which is connected to the lack of long-range order above 
T = in equilibrium correlations. Systems with textures 
in d > 2 have low-temperature ordered phases, and so are 
qualitatively different than ID and 2D systems. Indeed, 
we know of no scaling violations in systems with textures 
in d > 2. 



V. COMPARISON WITH PREVIOUS RESULTS 

In this section we review the current understanding of 
growth laws, and the evidence for our scaling assump- 
tions, for various phase-ordering systems. Other refer- 
ences can be found in reviews for scalar systems p],ppl| 
and in the review of Bray Q. The numerical simula- 
tions and experiments seem to be able to check whether 
a particular growth law (with or without logarithmic fac- 
tor) is reasonable, or to determine a growth-law exponent 
within about 10% without any theoretical input. The 
scaling of T(fc,t) is usually not considered. Within these 
constraints, existing results are consistent with the pre- 
dictions of the Energy-Scaling approach, as we discuss in 
more detail below. 



A. Non-conserved systems 

Scalar non-conserved systems, in agreement with our 
results, have growth laws L{t) ~ i^/'^ as long as the 
growth is driven by the curvature of domain walls, i.e. for 
spatial dimension d > 1 ||^,|l^ . This growth law has been 
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confirmed experimentally in two and three dimensions, 
see for example the work by Mason et at on twisted 
nematic liquid crystals ]T^ , and the work by Shannon 
et al. on Cu^Au [ p^ . In addition, simulations, such as 
Monte-Carlo (RG) studies also confirm this growth 
law. The scaling of the equal-time correlation function 
( pi] ) also has extensive confirmation, for example see the 
review by Furukawa In one dimension, scalar sys- 

tems exhibit growth driven by exponentially suppressed 
interactions between domain walls. This case has been 
treated approximately by Nagai and coworkers 0| who 
found a logarithmic growth law, and their solution and 
gowth law can be shown to be asymptotically correct 

Vector non-conserved systems have not been as well 
understood theoretically. We predict L ^ t^l"^ for all 
cases satisfying scaling, except for the 2D XY model 
where we expect L ~ (t/lnt)^/^. This has been con- 
firmed in the limit as the number of components n — > oo 
||9|,Q. The ID XY model has two time-dependent 
length scales A modified version of the Energy- 

Scaling treatment may still be used, however, and the 
growth laws determined if an independent relation be- 
tween the length scales can be established by physical 
arguments [E8|. Experimentally, liquid-crystal XY-like 
systems seem to exhibit L ^ t^/^ in two ||21|] and three 
pl| , p2t dimensions. 

In ID, simulations for systems without topological de- 
fects, with n = 3, 4, and 5, obtain a single length scale 
L ~ t^l"^ ^ while simulations for n = 2 are consistent with 
the scaling violations mentioned above 



In 2D, careful simulations of XY models were consis- 
tent with logarithmic factors, L ~ (t/lnt)^/^, in the 
growth law These factors manifest themselves as a 
systematic curvature in power-law fits of earlier numeri- 
cal studies 1^,^. However, tentative scaling violations 
in XY models have been found by Blundell and Bray 
p]| , when the length scale was extracted from the defect 
density rather than directly from the correlation func- 
tion. Similar inconsistencies are seen in simulations of 
uniaxial and biaxial 2D nematics (n = 2) by Zapotocky 
et al jS^]. No violations are seen when the correlations 
are collapsed with a length-scale extracted from the cor- 
relation function itself (see e.g. |^^). Measurement of 
r(k, t), the time-derivative correlations, may help to re- 
solve the scaling, or absence of scaling, of these systems. 
For n > 3 both hard spin {\(f\ = 1) and soft-spin 
p3[ simulations seem to find L t^^^. For n — 3 the 
hard spin simulations find L ~ t-^^^ from the energy den- 
sity (e L^^), but no scaling ||3^], which is confirmed 
by more extensive simulations of both hard and soft-spin 
texture systems [ ]69| . [Recent work by Toyoki p^ ] , which 
claimed scaling and L ^ t^/"^ in this system, is hampered 
by early-time transients.] 



pi[ . The suppressed growth law may be related to pin- 
ning effects due to the discrete dynamics — a possible 
explanation |3^,^^ for anomalously small growth expo- 
nents seen in simulations done at T = 0. Alternatively, 
corrections to scaling may be responsible. Simulations 
by Toyoki found L ~ f^/^ for 3D Heisenberg {n — 3) sys- 
tems . Scaling is confirmed by Blundell and Bray , 
though again they obtained growth exponents slightly 
less than 1/2. 

Our results derived from the scaling assumption are 
consistent with the bulk of the evidence. Further work 
needs to be done to measure T(k, t) in 2D XY systems to 
resolve the issue of scaling, and to check the deviations 
from t^/^ growth seen in d > 2 vector systems. 



B. Conserved systems 

We confirm that scalar conserved systems with fi — 2 
have the standard growth law oi L i^/"^ ]l7| , p8[ , pro- 
vided the spatial dimension d > 1 so that the evolution of 
domain walls is curvature driven. This result is also ob- 
tained by RG arguments [Q and is seen experimentally, 
for instance in the binary alloy Mno.eyC'uo.sa [ p9| , as well 

As in the non-conserved case. 



as in simulations |16 



0.45±0.01 



has 



In three dimensions, a growth law L 
been observed for the XY model |^,^,^,^ — less than 
the i^/^ prediction. However the system seems to scale 



scaling of ^(k, t) has been confirmed ||ll[ but the scaling 
of r(k, t) has not been considered. Kawakatsu and Mu- 
nakata |lO| have considered the ID case, in which there 
are no curvature effects, and found log arith mic growth in 
agreement with our results in section III C . Scalar mod- 
els with generalized conservation laws (p > 0) have also 
been considered, the so called "noninteger derivative" or 
"long-range exchange" models ||]. For scalar systems, 
growth laws in agreement with our results are found, 
both through computer simulations and by considering 
the evolution of small defect features. 

For vector conserved systems, with /i = 2, we predict 
growth laws of L ~ t^^^ — except for XY systems with 
d > 2 where we expect L (ilnf)^/*. This agrees with 
RG arguments, which, assuming standard scaling, found 
growth exponents of 1 /4 for all n > 2 ^] (we note again 
that the RG treatment does not determine logarithmic 
factors). Coniglio and Zannctti found scaling violations 
for n — oo | |59[ ], so that our treatment docs not apply to 
that case. However, Bray and Humayun | |73t show how 
standard scaling could be recovered for any finite n. 

Simulations of the conserved ID XY model find a scal- 
ing violation t hat i s consistent with the discussion at the 
end of section IV C 1 28] . Simulations of the 2D XY model 
|p6| obtain L ^ i /\~though small but systematic scal- 
ing violations seem to be indicated by the data. These 
"violations" could just be strong corrections to scaling, 
related to the corrections of order 0(1/ In Z) expected in 
the growth-law (from the cancelled logarithm from the 
energy-density in (p^). Recent simulations by Puri et 
al. found good scaling for both d = 2 and d — 3, 
with growth laws consistent with the forms L ^ t^/^ 
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(d = 2) and L ~ (tlnt)^/'^ (d = 3) predicted here. For 
the three-dimensional XY model, earlier work by Siegert 
and Rao |3^ also found good evidence for standard scal- 
ing. Although the growth law was originally interpreted 
as L ~ [Q, the data are much better fitted by the 
Energy-Scaling prediction L ^ (ilni)!/"* ||§. 

Again, our results are consistent with previous work, 
particularly the RG results for the power-law factors in 
the growth laws. Conserved 2D systems with textures 
(n = 3) should be investigated. In addition, the possible 
scaling violations in the conserved 2D XY system need 
to be resolved, perhaps by measuring ^(k, t). 



C. Long-range systems 

For systems with attractive long-range interactions 
( pl| ) , Bray Q] treated the conserved case with an RG ap- 
proach, and the non-conserved scalar case with physical 
arguments. Our results are in agreement, with additional 
logarithmic factors in marginal cases, for scalar and vec- 
tor systems for < cr < 2. Numerical work by Hayakawa 
et al. [|l2| for non-conserved 2D scalar systems roughly 
agree with our results for a = 0.5, 1.0, and 1.5, but more 
sensitivity is needed to test the logarithmic factor for 
a ~ 1. Analytical work by Hayakawa et al. found a 
breaking of scaling for conserved dynamics in the spheri- 
cal limit, similar to that seen in the short-range case 
Scaling violations in non- conserved ID scalar systems for 
cr < 1 have been proposed by Lee and Cardy EQ]. How- 
ever, no scaling violations in either S'(k, t) or T(k, t) are 
seen in more extensive simulations of ID scalar systems 
[^, which find growth laws in agreement with our pre- 
dictions for cr = 0.5, 1.0, and 1.5. 



VI. DISCUSSION 

Our approach does not depend on the initial condi- 
tions of the system If scaling is obeyed, critical 

(^<?^ = 0) and off-critical (^^^ 7^ 0) quenches will have 
identical growth laws. For non-conserved systems, off- 
critical quenches are found to break scaling [^,0 , since 
the order-parameter saturates, and so our growth laws do 
not apply. Indeed, exponentially growing length-scales 
are seen Conserved systems, on the other hand, 

have the same growth law for critical and off critical 
quenches, in agreement with earlier predictions Jl8| , ^ . 
This has been numerically confirmed in the scalar case 
jTSf . Initial conditions with long-range correlations (cor- 
responding to, e.g., quenches from a critical point) will, 
similarly, have the same growth laws if scaling holds . 
In general, various classes of initial conditions may affect 
the existence of scaling and the form of the scaling func- 
tions (|l^) and (^3|), but will not change the growth laws 
if scaling holds. Of course, the appropriate defect struc- 
ture must be applied. For example, if the 2D XY model 



is quenched from below Tkt then the bound vortex pairs 
quickly annihilate and the scaling configurations have no 
topological defects (for a quench to T = 0). In this case 
scaling is obeyed and the growth law is. L ^ t^/'^ [Q, 
as expected for a defect-free quench, without the loga- 
rithmic factor of 2D XY systems quenched from above 
Tkt- 

We derive our results at zero-temperature but we ex- 
pect our growth laws to describe quenches to all temper- 
atures below Tc- We do not address quenches in which 
thermal noise is essential, such as systems with static dis- 
order [|o| , or quenches to a T > critical point |^l| . For 
systems with long-range order only at T = 0, such as 
XY and scalar systems in one-dimension, our approach 
applies only at T = 0. For small T, such that the equi- 
librium correlation length ^(T) of the disordered phase 
is large, we expect our growth laws to apply as long as 
L{t) < f(r). The 2D XY model quenched from high- 
temperature to < T < TxT develops power-law order 
and invites further study because the low-temperature 
behavior is not described by a single T = fixed point 
but by a line of fixed points. Numerical work for this case 
by Yurke et al. at < T < Tkt, found a growth law 
in agreement with our T = result, suggesting a similar 
growth-law for all temperatures below Tkt- 

Our treatment does not depend on the details of the 
potential V{(j)) in the energy- functional. It is only the 
symmetry properties of the potential's ground-state man- 
ifold that determine the growth law, since the ground- 
state manifold determines the defect structure, which in 
turn determines the asymptotic structure at kL ^ 1 
through simple scaling arguments. The Energy-Scaling 
approach can be applied to systems with more compli- 
cated order parameters than n-component vectors. All 
we need is the existence of some short or long-ranged 
"elastic" energy (cr), a conservation law (/x), and the de- 
fect structure, if any. The effective n corresponding to 
a defect type is the one that determines the scaling of 
its core volume density Pcorc ^ i^". The defect type 
with the smallest n, and hence greatest core volume, 
will determine the growth law, since it will provide the 
dominant contribution to any UV divergences in the en- 
ergy density or dissipation integrals. Systems with non- 
abelian symmetries, such as cholesteric liquid crystals, 
are treated in the same way, since our approach is predi- 
cated on the energetics of the system rather than on the 
detailed nature of the dynamics. This assumes that scal- 
ing holds, so that all defect types scale with the same 
growth law. 

For example, in bulk uniaxial or biaxial nematic liq- 
uid crystals the existence of string defects, with Pstiing ^ 
L~^, leads to n = 2. Using this in equations (Il3) and 
( [TtI ) with no conservation law implies a L ^ t^/^growth 
law. This is consistent with recent experiments [[2l[|22| 
and simulations |30|. We neglect the point defects, since 
they do not dominate the asymptotics or the energetics. 

Similarly, in Potts models the existence of domain walls 
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leads to the same growth laws as for scalar (n = 1) sys- 
tems, consistent with the L ~ i^^^ growth seen in non- 
conserved systems [Q, and the L ~ t^/^ growth sug- 
gested by numerical studies of conserved systems [q3|. 
The existence of vertices in Potts models, or e.g. in 
clock models , will not change the form of the asymp- 
totic growth law. At late times the energy density and 
dissipation will be dominated by the change of domain 
wall volume, with pwaii ^ L^^, rather than by the ver- 
tices, with Pvcrticcs ~ L^^ . In general, the prefactor or 
timescale of the growth law will depend on the details of 
the model, but the exponent and any logarithmic factors 
of the growth law will be universal. 

Our approach is restricted to systems which are gov- 
erned by dissipative dynamics and a simple energy- 
functional of the form (|l|) or (pT|). It would be inter- 
esting to speculate on the growth laws in the late stages 
of defect elimination in patterned systems, with compet- 
ing short and long-range interactions. However, while 
the defect types of the patterned structures (disclinations 
etc..) can be clearly identified, it is not at all clear what 
conservation laws or long-range forces apply to the effec- 
tive order parameter of the patterned structure, or even 
whether temperature is relevant to phase-ordering [ p5| . 
This is a promising direction for further research. 

For the systems where we know the "collapse" laws 
lit), of the size of an isolated defect for a given time to 
collapse, they are the same form as the scaling growth law 
L(t), of the phase-ordering length-scale for a given time 
after the quench. This agrees with our naive expectation 
for scaling systems that the form of the collapse law will 
be the same as the growth law — the intuitive picture 
is of collapsing defect features leaving "voids" which set 
the growing length-scale, L{t). 

We can say little about the growth laws in systems that 
break our scaling assumptions, apart from our discussion 
about the dependence at the end of section |l|. Some 
systems can be explicitly shown to break scaling, such 
as ID XY models or conserved spherical systems 
Are there other systems that break scaling? We 
cannot answer this by examining only the equal-time cor- 
relation function, since the scaling of the latter does not 
imply the scaling of T(k, t). There are two paths to take. 
The first is to explicitly find scaling violations in a sys- 
tem. A possible example is given by the equal-time corre- 
lations in the 2D XY model, both for the non-conserved 
pl| , ^ and for the conserved (apparent in figure 3 of |26|) 
cases. The second method is to find a growth law that is 
in striking disagreement with the our predictions, such as 
in the 2D n = 3 system . Since our approach is based 
on the scaling assumption, any disagreement implies a 
scaling violation. Conversely, to demonstrate scaling the 
correlations must be measured directly since agreement 
between the observed and predicted growth laws is nec- 
essary but not sufficient to demonstrate scaling. 



VII. SUMMARY 

In summary, by focusing on the energetics, rather than 
the detailed dynamics of the system, we obtain growth 
laws for phase ordering ( summarized in figures 

and |, 

as well as table IV). This leads to a powerful, unified, 
and physical approach to determining growth laws that 
rests solely on the existence of scaling. We call this the 
"Energy-Scaling" approach. It can be used for any sys- 
tem with purely dissipative dynamics. Any disagreement 
of a growth law from our predictions indicates the break- 
ing of scaling of at least one of the two-point correlation 
functions. In particular, we stress the importance of the 
time-derivative correlation function. 

In addition to the growth laws, we determine the col- 
lapse law for isolated non-conserved textures {n = d-\-l) 
in more than two-dimensions, and use that to determine 
the generalized Porod's law (Q) and asymptotics for 
r(k, t) (^9|) — these follow novel power-laws due to anni- 
hilating textures. We also treat the dynamics of isolated 
defects in scalar systems and 2D XY systems. When scal- 
ing holds, we find that the growth law matches the form 
set by collapsing isolated defect features, of the charac- 
teristic scale L, within a larger phase-ordering system. 
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FIG. 2. The growth law of the length scale L{t) for attrac- 
tive long-range interactions with a fixed sigma, < a < 2 for 
systems which satisfy our scaling assumptions. The growth 
laws are shown for various number of components n and con- 
servation laws fi. For n = d = 1 < a < 2, L{t) ~ ti/(i+<^) 
[^. For a = 2 and n = d < 2, see table Scaling sys- 
tems with no topological defects (n > 1), or with textures 
(n = d + 1), will have L ~ 

FIG. 3. The construction of a d -I- 1 component texture 
in a d-dimensional system by stereographically projecting a 
d-sphere onto the system (shown is a plane containing x 
and the two-poles p± of the d-sphere). The order-parameter 
(/f>(x) = R, where the radial vector R meets the d-sphere at 
the intersection with a line joining x with the pole p+ of the 
sphere. The scale of the constructed texture, X, is twice the 
radius of the d-sphere. 

TABLE I. Results of the energy dissipation and the domain 
size dynamics of an isolated scalar domain of scale I <Si L with 
d > 1. The time left before domain annihilation is t. 



TABLE II. Evolution of a pair of defects of separation 
/ < L for 71 = d = 2. 

TABLE HI. The average square- velocity, (i"^) , of defect 
cores which contribute to T(k,t) at kL 3> 1 for scalar and 2D 
XY systems. The dots indicate time derivatives. 

TABLE IV. Growth laws, L{t), for short-range systems 
with n — d < 2. 



FIG. 1. The growth law L(t) of the length scale for 
phase-ordering systems which scale. The growth laws are 
shown for various number of components n, where n < d, 
and conservation laws, jj.. The black circles correspond to 
possible physical systems. For n = d = 1 and n = d = 2 
the growth-laws in table ^ apply- Scaling systems with no 
topological defects (n > d + 1), or with textures (n = d -I- 1), 
wiU have L 
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